Tunneling Spectroscopy and Vortex Imaging in Boron-Doped Diamond 
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We present the first scanning tunneling spectroscopy study of single-crystalline boron doped 
diamond. The measurements were performed below 100 mK with a low temperature scanning 
tunneling microscope. The tunneling density of states displays a clear superconducting gap. The 
temperature evolution of the order parameter follows the weak coupling BCS law with A(0) //csT c ~ 
1.74. Vortex imaging at low magnetic field also reveals localized states inside the vortex core that 
are unexpected for such a dirty superconductor. 

PACS numbers: 74.50.+r, 74.25.Qt, 74.78.Db 



The recent discovery of a superconducting transition in 
boron-doped diamond provides a new interesting system 
for the study of superconductivity in doped semiconduc- 
tors 0. In the broad band gap of diamond (~ 5.5 eV), 
boron impurities introduce an acceptor level with a hole 
binding energy of ~ 0.37 eV and lead to a metallic state 
above a critical boron concentration in the range of a 
few atoms per thousand. According to the doping level, 
the superconducting transition is observed with a critical 
temperature varying between IK and 10 K 0,0. 

Because of the strong bonding states and the hole dop- 
ing, various theoretical studies have stressed the simi- 
larity between diamond and the superconductor MgB^ 
jj Li U] • The origin of the remarkable superconductivity 
in MgF>2 at T c = 39K is now well understood to be a 
result of phonons coupled to holes in the two dimensional 
a-bonding valence bands. Although the band structure 
of diamond is three dimensional, the superconductivity 
in boron doped diamond could also be due to the cou- 
pling of a few holes at the top of the valence band to the 
optical bond-stretching zone-center phonons 00 as well 
as to boron-related modes 0, . Numerical calculations 
of the electron-phonon coupling by the density-functional 
supercell method are claimed to yield a calculated T c in 
good agreement with experiments for at least one boron 
concentration 0, thus supporting that boron-doped di- 
amond is a phonon mediated superconductor. 

Another theoretical approach proposed by Baskaran 
considers diamond, a broad band insulator, to be an 
appropriate vacuum state for the boron subsystem; elec- 
tron correlations then drive a Anderson-Mott insulator 
to resonating valence bond superconductor transition . 
Because of the proximity to the boron critical density 
of the insulator-metal transition, the superconductivity 
takes place within the impurity band. In the disordered 
lattice of boron acceptors, such a strongly correlated im- 
purity band is expected to give rise to an extended s-wave 
superconducting state. 

In this letter we report on low temperature tunnel- 
ing spectroscopy and vortex images of superconducting 



hole-doped diamond films. These measurements which 
probe the quasiparticule excitations near the Fermi en- 
ergy, show a clear BCS Local Density of States (LDOS) 
consistent with weak coupling. Contrary to what is ex- 
pected for a dirty superconductor [lQ, [u| , a significant 
density of localized resonant states is found below the 
gap in the vortex core. 

The sample studied was a selected boron-doped single- 
crystalline diamond epilayer synthesized at 820°C by mi- 
crowave plasma-assisted decomposition (MPCVD) of a 
gas mixture composed of Hydrogen, Methane and Dib- 
orane with a controlled boron to carbon concentration 
ratio. A 001-oriented type lb diamond substrate was ex- 
posed to a pure hydrogen plasma prior to diborane-free 
growth of a 500 nm-thick buffer layer. Optical and trans- 
port measurements showed that these steps were critical 
for achieving a high quality boron-doped diamond epi- 
layer. A boron atomic density varying between n# = 



3.1 xlO 21 



at the surface and = 2.1 x 10 21 



at the interface, as well as an epilayer thickness of 75 nm 
were derived from a quantitative comparison to a Sec- 
ondary Ion Mass Spectroscopy (SIMS) profile measured 
in a B-implanted diamond crystal with a known peak 
boron concentration. Because the density and effective 
masses of the charge carriers are not clearly established 
0, transport caracteristics such as the Fermi velocity 
Vf and mean free path I, can only be roughly estimated. 
For a carrier density deduced from the boron density 
and a sheet resistance of 83f2/n at 4.2 K, I is estimated 
within the free electron model to be less than 1.3 nm. The 
surface topography obtained with the scanning tunneling 
microscope (STM) is displayed in Fig. [I] The parallel 
600 nm wide stripes reflect the vicinal surface structure 
of the substrate and can be observed everywhere with the 
same orientation. The rms micro-roughness is 1.8 nm and 
the in plane characteristic size of the surface corrugation 
is about 75 ± 25 nm. 

Spectroscopic measurements were performed with a 
home-built Scanning Tunneling Microscope (STM)cooled 
down to a base temperature of 50 mK in a dilution re- 
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FIG. 1: Topography image (1.5 x 1.5 /im 2 ) of the sample 
obtained at 70 mK. The rms rugosity is 1.8 nm. 



frigerator. The sample was scanned with a Pt-Ir cut tip. 
In order to probe the LDOS, a small AC modulation of 
10 jiV rms was added to the sample-tip DC bias volt- 
age V and the differential conductance G(V) = ^(V) 
obtained with a lock- in amplifier. At any position above 
the sample, V could be ramped and the resulting G(V) 
curve giving the LDOS around the Fermi energy mea- 
sured. In addition, a static magnetic field perpendicular 
to the sample surface could be applied with a supercon- 
ducting coil surrounding the STM. 

Tunneling spectroscopy performed at different loca- 
tions revealed the surface to be superconducting with 
very little spatial inhomogeneity. In Fig. Owe show two 
representative differential conductance curves performed 
at 70 mK with a tunnel resistance of 20 MQ and without 
any averaging. Most of the experimental spectra can be 
well reproduced by a theoretical BCS density of states as 
can be seen in Fig. [21(a). The fit gives a superconducting 
order parameter A = ( 285 ± 2 ) fieV. A thermal broad- 
ening with an effective temperature T e ff = 235 mK in 
the LDOS calculation indicates that the energy resolu- 
tion of our STM is probably limited by unfiltered elec- 
tromagnetic radiation which heats the electron bath. No 
additionnal parameters to describe pair breaking were 
needed. Although the spectrum was almost identical 
everywhere on the scanned surface, we also obtained 
smaller values of A in several independent experiments 
on this sample (for example, see Fig. 0(b)). This slight 
dispersion of the gap value can be either the result of spa- 
tial macroscopic inhomogeneity of the doping or be the 
consequence of variation in the chemical cleaning of the 
surface prior to each run. The temperature dependence 
of the order parameter A(T) displayed in Fig. [2| is well 
described by the BCS theory with a critical temperature 
of 1.85 K very close to the superconducting transition 
temperature of 1.9 K obtained by ac susceptibility and 
transport measurements. This clearly demonstrates that 
boron-doped diamond is well described by s-wave BCS 
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FIG. 2: Top : experimental normalized tunneling conduc- 
tance measured at 70 mK (open circle). Solid lines corre- 
spond to BCS fits with : (a) A = 285 fieV and T e // = 
235 mK; (b) A = 260 fieV and T eff = 175 mK. Bottom : 
(c) temperature dependance of the BCS gap (open squares) 
compared with a BCS law with T c = 1.85 K (solid line). 



superconductivity with a measured ratio ~ 1.74 

characteristic of weak-coupling hjj , as expected from 
theoretical calculations jj, la, Lij • Nevertheless we also 
found a small number of spectra as shown in Fig (b) 
which present broader gap edges. In order to explain 
them, further investigations of the LDOS for different 
crystallographic orientations are needed to investigate a 
possible anisotropy of the order parameter. 

In a magnetic field, for a type II superconductor a 
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FIG. 3: Top : vortex images (1.5 x 1.5 /im 2 ) at two different 
magnetic fields : (a) 1200 Oe, (b) 1900 Oe. Bottom : two di- 
mensional autocorrelation function calculated from the above 
vortex images. 



mixed state forms where magnetic vortices each carry- 
ing one flux quantum </>o penetrate the sample. In the 
absence of pinning and for an isotropic order parameter, 
the repulsive interaction between flux lines give rise to 
a periodic lattice of vortices with a triangular symmetry 
: the Abrikosov vortex lattice 14]. A vortex represents 
a singularity with the superconducting order parameter 
suppressed at its center. In the clean limit, i.e. £o/l <C 1 
(£o = Ja(o) ^ s tne coherence length) confinement of 
quasiparticles inside the vortex yields discretized energy 
levels separated in energy by A 2 /Ep 15]. Such bound 
states were first observed in NbSe2 as a Zero Bias Con- 
ductance (ZBC) peak [l^. STM studies of vortices in 
high temperature superconductors have revealed no ZBC 
peak but additional quasi-particles states at constant en- 
ergies above and below the Fermi level 0, ^| . In the 
less theoretically studied dirty limit case, i.e. ^o/l ^> 1, 
a mixing of the core bound states by disorder results 
in a flat normal-metal-like LDOS predicted ^(j and ex- 
perimentally observed in Nb\- x Ta x Se2 at 1.3 K [l]| . 
The temperature independent vortex size is then equal 
to the dirty coherence length £ s = (hD /27t/c j bT' c ) 1//2 , 
where D is the diffusion coefficient. For our sample, up- 
per critical field measurements yield a Ginzburg-Landau 
coherence length ^ql = 15 nm from the expression 
H C 2 = 0o/(27r£(? L ) .2]. With a mean free path Z ~ 1.3 nm 
and a BCS coherence length £o ~ 240 nm evaluated from 
the relation £gl = 0.855^ol, boron-doped diamond is 
clearly in the dirty limit with £ s = (2/tt)( i gl — 9.5 nm. 



We have imaged the vortex arrangement by combining 
topography and spectroscopy measurements in a mag- 
netic field applied perpendicular to the diamond film. 
The sample-tip DC bias was set to a voltage V close 
to the BCS superconducting coherence peak (eV = A). 
Then, the presence of vortex is revealed by a decrease of 
jjy which gave a dark spot in the differential conductance 
image. In Fig. [31(a) and (b) we show two 1.5 x 1.5 jam 2 
vortex images obtained at about 100 mK and induced 
by a field H of 1200 Oe and 1900 Oe respectively. The 
vortex positions were very stable and insensitive to the 
scanning. In order to reduce the noise, these images were 
convoluted with a disc of diameter charateristic of the 
lentgh scale over which the contrast changes occur in the 
raw image. Its radius R is not a direct measurement of 
the coherence length £ s because the decrease of the coher- 
ence peak begins well before crossing the vortex core [l0| . 
Indeed, we measured R = 40 ± 5 nm and R = 30 ± 5 nm 
for a magnetic field of 1200 and 1900 Oe respectively. 
This dependency can be explained by a field enhanced 
vortex- vortex interaction (l0( . 

Contrary to the perfect Abrikosov vortex lattice, the 
spatial distribution of vortices is strongly disordered even 
though a local hexagonal arrangement persists. In order 
to characterize it, the autocorrelation function defined 
by G(5x, Sy) = J I(x + 5x, y + 5y)I(x, y) dx dy was calcu- 
lated. Here, I(x, y) is the local degree of brightness. G(S) 
gives the probability of finding a similar brightness in the 
image for a spatial shift equals to 5 = Sy) from any 
point. The results are shown on Fig. E|(c) and (d) for the 
vortex images Fig. 03 (a) and (b) respectively. Although 
the lattice appears strongly disordered, the sixfold sym- 
metry of first rings of neighbors in the autocorrelation 
pictures indicates a persistent orientational order for any 
site of the lattice. The mean distance between vortices is 
equal to 110 ± 5 nm and 140 ± 1 nm, in good agreement 
with d/\ = j^-) 1 ^ 2 expected for a triangular lattice. 
Moreover, for the higher field we can see a rotation and 
a slight deformation of the ring of first neighbors accom- 
panied by the onset of a diffuse striation in a direction 
which tends to align with the vicinal stripes of the sur- 
face visible on Fig. ^ This alignment is consistent with 
a stronger pinning by these stripes at higher flux density. 

In order to map the LDOS inside a vortex we ac- 
quired spectra every 0.7 nm along lines equally spaced 
by 2.2 nm. The exact localization of the core is a deli- 
cate problem in a dirty superconductor since disorder is 
expected to smear out the ZBC peak by scattering the 
localized quasiparticule states and consequently yield a 
flat normal- metal-like LDOS 0,^1- Here, instead, we 
surprisingly obtained numerous localized resonances at 
non zero energies inside the gap while crossing a vor- 
tex along different lines. We consider a vortex core to 
be defined as the region containing the highest concen- 
tration of resonance peaks within our spatial resolution. 
The corresponding line of spectra is displayed in Fig. 
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FIG. 4: Tunneling spectroscopy at T — 53 mK through a 
vortex core induced by a magnetic field of 1800 Oe. Top : 
tunneling spectra acquired every 0.7 nm along a line crossing 
the vortex core which extends between the two arrows. A slid- 
ing averaging of two consecutive spectra was applied in order 
to attenuate high peaks for clarity. Bottom : two-dimensional 
view of the same data restricted to the vortex core and with- 
out averaging. The color scale as well as the contour lines give 
the normalized LDOS. Black arrows indicate the positions of 
the localized resonances. 



0] (top). Such a determination gives a vortex core ra- 
dius of about 10 nm consistent with £ s . More striking is 
the periodic-like spatial occurrence of these resonances as 
pointed by the black arrows in the two-dimensionnal view 
of Fig. 21 It is worth noticing that a periodic pattern of 
quasi-particles states has also been observed around vor- 
tex cores in Bi^Sr^ C aCu<2 Op+fi , another superconductor 
close to a Mott transition |l£j. The distance of about 
2.2 nm between each localized resonance is almost con- 
stant and present on most scanned line. Such a spatial 
modulation appears in the theoretical clean limit [2Q( at 
low temperature T/T c <C 1 as Friedel-like oscillations 
with a period of the order of the Fermi wave-length. 
Although the distance between each localized resonance 



is consistent with an estimation of \p < 1.6 nm in a 
free-electron model, our sample is not in the clean limit. 
Moreover, because of thermal smearing, quantum states 
at a energy of A 2 /Ep, would appear as a ZBC peak. On 
the other hand, these puzzling resonances could be asso- 
ciated to low-lying vortex bound states shifted away from 
the Fermi level by the impurity potential 0,[2^. Then, 
the spatial position and the energy of these resonances 
would be specific to the local configuration of disorder. 

In conclusion, we performed the first spectroscopic 
tunneling study and STM images of the vortex lattice 
in superconducting boron-doped diamond. Our results 
provide a local density of states in excellent agreement 
with the BCS s-wave theory with a ratio A//cbT c c± 1.74 
characteristic of a weak coupling superconductor. Under 
magnetic field, vortices which are arranged in a strongly 
disordered triangular lattice present unexpected reso- 
nances inside the BCS gap at non zero energies. Fur- 
ther investigations are needed in order to understand 
whether these magnetic field-induced resonances are a 
generic feature of any superconducting state neighboring 
a Anderson-Mott insulator. 

The authors are indebted to C. Cytermann (Technion 
Haifa) for measuring the SIMS profile of the sample anal- 
ysed in this work. 
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